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The differences between spin relaxation in grapfiene and in other materials are discussed. For 
relaxation by scattering processes, the Elliot-Yafet mechanism, the relation between the spin and 
the momentum scattering times acquires a dependence on the carrier density, which is independent 
of the scattering mechanism and the relation between mobility and carrier concentration. This 
dependence puts severe restrictions on the origin of the spin relaxation in graphene. The density 
dependence of the spin relaxation allows us to distinguish between ordinary impurities and defects 
which modify locally the spin-orbit interaction. 
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Introduction. Graphene is considered as a potential 
material for spintronics devices due to the weak spin- 
orbit (SO) interaction^ and long spin lifetimes^. One of 
the most intriguing features of spintronics in graphene 
is that the observed spin relaxation time is significantly 
shorter than the values estimated theoretically. A num- 
ber of recent experiments^^— investigate which spin re- 
laxation mechanism plays the major role both in single 
layer and in bilayer graphene. The D'yakonov-Perel^ and 
the Elliot- Yafet^i^ mechanisms have been discussed in the 
context of graphene^iiS. Experiments suggest that the 
main source of spin relaxation in single layer graphene 
is extrinsic, lending support to the Elliot-Yafet mecha- 
nism. Longer spin lifetimes have been reported in bilayer 
graphene than in single layer—, when the SO coupling in 
the bilayer is expected to be somewhat stronger—. 

The Elliot-Yafet mechanism takes into account the 
change in the spin polarization of a Bloch electron due to 
scattering by impurities, lattice defects or phonons. The 
Elliot relation establishes a linear relation between the 
spin relaxation time and momentum scattering time: 



a 



(1) 



where a can be interpreted as the spin-flip probability 
during a momentum relaxation event. Elliot deduced this 
relation by using a pcrturbative approach. Due to the SO 
coupling, Bloch states with well-defined spin polarization 
are not longer eigenstates of the complete Hamiltonian. 
In the case of conventional metals with a center of sym- 
metry, two degenerate states can be defined for each value 
of the momentum^: 



[ak(r)|t>+&k(r)|;)]e^' 
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(2) 
(3) 



where the coefficients a, b are lattice-periodic due to the 
discrete translation symmetry. These two states are con- 
nected by spatial inversion and time reversal symmetries 
and form a Kramers' doublet. Typically these states 
can be identified with spin-up and spin-down states be- 
cause typically |6| <C 1. Since the SO interaction cou- 
ples electronic states with opposite spin projections in 



different bands (in the case of graphene the SO interac- 
tion couples TT and a bands), perturbation theory gives 
\b\ fa Aso/^E, where AE is the energy difference be- 
tween the two bands involved. Usually, Aso ^ ^E, as 
in the case of graphene. 

The spin flip amplitude during the scattering by an 
obstacle with no spin degrees of freedom itself can be 
computed using the Born approximation, leading to a « 
<^|&P) where the symbol () expresses an average over the 
Fermi surface. These arguments are quite general and 
do not depend on the nature of the scatterers. Realistic 
calculations can be done in some cases, for instance in 
the case of III-V semiconductor a^^'^'^ . 

The relation ([T]) holds experimentally for most conven- 
tional metalsii^. As we discuss below, doped graphene 
is not an exception. However, unlike ordinary metals, 
the nature of the effective SO coupling acting on the 
graphene tt electrons, which are the relevant ones in 
what concerns to transport properties, makes the ra- 
tio Ts/Tp to depend strongly on the number of carriers 
through the Fermi energy. This result holds for differ- 
ent kind of defects, as it discussed next. A wide va- 
riety of experimentsi^ suggest a linear scaling between 
Tp and Ts, with independence of the carrier concentra- 
tion. Our analysis shows that this behavior cannot 
attributed to the Elliot-Yafet mechanism, opening the 
door to other extrinsically induced spin relaxation mech- 
anisms, such as a defects which modify locally the spin- 
orbit interactioni^iii. 

The model. In graphene, the intrinsic SO coupling can 
be neglected in comparison to the Rashba-like coupling, 
generated by perturbations which break spatial inversion, 
such as electric fields and ripples. If the perturbation 
changes slowly over scales larger than the lattice spac- 
ing, we can neglect intervalley hybridization^. Then, 
the Hamiltonian of the problem reads: 



A 
H = —ikvpo • V + — ((J X s)^ 



(4) 



The Rashba-like term breaks the spatial inversion sym- 
metry, and two degenerate eigenstates cannot be defined 
for a given momentum k. The Rashba-like term entan- 
gles spin and valley degrees of freedom, complicating the 



definition of the amount of spin relaxation in a scattering 
event. The Bloch eigenstates of (|3]) read: 
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where ^k 



arctan (fcy/fca;) and e^ . 



±^ 



(?it;i?|k|)2 + (-y) , where e denotes electrons. A simi- 
lar expression can be defined for holes by changing the 
sign of the second term. In what follows, we restrict the 
discussion to electrons. As we see, a spin direction can- 
not be uniquely defined for all momenta. When we take 
A = 0, eigenstates ([S]) are Bloch states with well-defined 
projection of spin over the direction of motion, that is, 
helicity ±. This is not strictly true when A 7^ 0, but 
in the spirit of the above Elliot's approach, we can iden- 
tify each of these eigenstates with chiral states ±. This 
is justified from the point of view of perturbation the- 
ory, since for carrier concentrations of interest we have 
A/ep ^ 1- Thus, the effect of the Rashba-like coupling 
can be interpreted as the energy splitting {'^ A) of bands 
with opposite chirality. 

Let's consider now scattering by a potential U (r) di- 
agonal in sublattice and spin degrees of freedom in the 
Born approximation. We study scattering in the chiral 
channels discussed above instead of the spin-up and spin- 
down channels as in the case of the Elliot's approach. 
This restriction complicates the definition of the amount 
of spin relaxation. To illustrate this, it is useful to cal- 
culate the scattering amplitudes in these channels in the 
absence of SO. Assuming and incoming Bloch state with 
energy e = hvpk and positive chirality, it is easy to see 
that in that case (see Supplementary Information): 



(6) 



where C/q is the Fourier transformation of the scatter- 
ing potential evaluated at the transferred momentum 
k' — k, and 6 is the angle between the outcoming k' 
and incoming k momentum (see Fig. [T]). If we repeat 
the calculation in the spin-up and spin-down channels 
assuming an incoming state with spin up, then we ob- 
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tain /t = - (hvp) y-i^Uq {l + e *'') and /j, = 0, since 
in the absence of SO there is no spin-fiip. The scat- 
tering amplitude /_ is not zero in general (except for 
forward scattering) so that it cannot be related with a 
cross-section for a spin-fiip process. 

In order to study scattering in the chiral channels de- 
fined by the Rashba coupling we define the probability 
for a spin-flip process from the changes in the scattering 
in both chiral channels due to the presence of the SO cou- 
pling. We follow the approach of Ref. [l| adapted to the 





FIG. 1: a) Sketch of scattering by a potential U (r) in the 
chiral channels defined in the text, b) Sketch of scattering by 
a boundary. 



calculation within the Born approximation. We define 
the quantity: 



3(0) 
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where the superscript A (0) indicates the presence (ab- 
sence) of the Rashba-like coupling. This quantity van- 
ishes when A = 0, and it can be interpreted as a mea- 
sure of the amount of spin relaxed in the direction defined 
by 6*. As in the case of the Elliot's approach, the total 
amount of spin relaxation during a scattering event can 
be defined as the average of this quantity over the Fermi 
surface: 



S={S{6))^^jd6S{e,e = eF) 



(8) 



In Ref. m the relation S ~ A/ep was deduced in the 
case of weak scatterers. Here we show that this rela- 
tion is general, and it does not depend on the nature of 
the scatterer, including strong scatterers or other impu- 
rity potentials which cannot be treated in the Born ap- 
proximation, where the value of S cannot obtained from 
perturbation theory. The Born approximation suffices, 
however, to show how this behavior is implied by the na- 
ture of the SO coupling in graphene with independence 
of the precise scattering mechanism. It is not difficult 
to compute exactly f^ [0) (see Supplementary Informa- 
tion), but the picture provided by perturbation theory is 
enough to illustrate this behavior in the doped regime. 
For an arbitrary scatterer, the value amplitude f±{d), 
as defined in (|6]), requires the use of non perturbative 
methods. The difference f^{0) — f±{6), however, can be 
obtained by expanding in powers of A/ep. This is easy 
to see in the Born approximation, where the substitution 
e — )- e± A/2 in expressions © has to be made in order to 
obtain f±{9). An expansion in powers of A/e is well de- 
fined, and it implies that S (0) ~ A/e, independently of 
the scattering potential, [/q, which factorizes in expres- 
sions ^. Assuming this behavior, the Elliot relation for 
graphene can be easily found. After Ncoi collisions, the 



change of spin polarization is of the order of ^/NcoiS. 
Dephasing takes place after a time Ts = NcoiTp, when 
\/NcoiS ^ 1. Hence we obtain the relation: 



,2 'P 



(9) 



This is the Elliot relation for graphene. As one can 
see, the ratio Ts/Tp depends on the carrier concentration 
through the Fermi energy. In what follows we compute 
exactly the amount of spin relaxation S for different kind 
of scatterers, generalizing the relation ([9]). 

Results for different scatterers. The scattering ampli- 
tudes ([6]) can be calculated exactly in the presence of the 
Rashba-like coupling (sec Supplementary Information). 
In the case of weak scatterers, we consider as scattering 
center a isotropic potential C/ (r) = Vd{r — R), where 
i9 (r — i?) is a step function. In the case of Coulomb scat- 
terers the scattering potential reads U (r) = —hvpajr. 
Note that |a| < 1/2, in other case the solutions of the 
Coulomb problem oscillate very fast and have no well- 
defined limit as r — > 0, which corresponds to the Dirac 
vacuum breakdown (the continuum description in terms 
of the Dirac Hamiltonian is not valid)^^. The results are 
shown in Fig. [51 In both cases S ^ ^/ep- 



In order to study spin relaxation during scattering by 
a boundary we have to adapt the definition of S. We 
are going to consider a zig-zag termination for simplicity, 
since it defines the most general boundary conditional^. 
We consider as incoming wave a Bloch state ^k.+ with 
energy e, forming an angle (j) with the direction perpen- 
dicular to the boundary. As it is deduced from Fig. [TJ 
7r/2 -\- a ^ (j), where a ~ arctan {ky / k^) ■ Two outgo- 
ing Bloch states exist satisfying conservation of energy 
and momentum in the direction parallel to the boundary. 
Then, the outgoing wave can be written as the superposi- 
tion ^ojit = ''i^k+,-1- +?'2^k_,-, where k_|_ (k_) forms an 
angle a' {cr") with the direction defined by the boundary 
(see Fig.[T]), and |k±| = fc± = {hvpY Ve^ =F eA. We can 
define the amount of spin relaxed in the direction defined 
by as: 
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where the superscript refers to the refiection coefficients 
in the absence of the SO coupling. As before, the amount 
of spin relaxed by the boundary can be defined as the av- 
erage, S = {S (0)) = i /^^/2 d(j)S ((/), e = ef)- By impos- 
ing zig-zag boundary conditions we obtain the following 
expressions for the refiection coefficients: 
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FIG. 2: S as a function of the carrier concentration for A = 1 
nieV (in blue) and A = 0.5 meV (in red), a) Weak scatterers 
(i? = 1 A and Vo = 0.1 eV). b) Coulomb scatterers. Tliese 
results are obtained by computing exactly the scattering am- 
plitudes f± and evaluating numerically Eq. ((HJ. Insets: S for 
A = 0.1 meV and kp — (hvp)' ep — 0.01 A plotted as func- 
tions of kp^ and A respectively. A clear linear dependence is 
showed, as it is argued in the text. 



Besides this, from conservation constrains we have a' = 
—a and cos a" ~ j-^ cos a. To first order in the SO cou- 

pling, a" = -cr - f cot cr + ©(A^/e^). Then, to first 
order in the SO coupling, the refiection coefficients in 
terms of the angle (j) read: 
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(13) 
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The amount of spin relaxation can be estimated as 5 = 
-2^ . This expression fits the exact result rather well, see 

Film 

In the case of strong scatterers, such as vacancies, the 
Born approximation fails. As in the case of boundaries, 
we need to extend the definition of 5*. Strong scatterers 
can be described as a circular void of radius of the order 
of the lattice constant. We can exploit the cylindrical 
symmetry of the problem by using the decomposition of 
the eigenstates of (|4]) into partial waves with well-defined 
generalized total angular momentum J = lz+(7z/'2+Sz/2, 
which is actually a global symmetry of the problem, 
where l^ is the third component of the orbital angular 
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FIG. 3: S as a function of the carrier concentration for A = 1 
meV (in blue) and A — 0.5 meV (in red). The dots correspond 
to the numerical evaluation of S, the continuum line to the 
analytical estimates of the text, a) Scattering by boundaries, 
b) Strong scatterers {R = 1.4 A). 



momentum operator Iz — —i{xdy — ydx)- For each in- 
coming cylindrical wave with energy e, there are two re- 
flected waves with the same energy. A quantity analo- 
gous to Eq. ([TU|) can be defined, considering scattering 
in all channels with well-defined J (see Supplementary 
Information). The amplitude for spin scattering can be 
calculated analytically 5* w |^ which fits very well the 
numerical evaluation of S, see Fig. |3l 

The case of clusters of impuritics^i is studied within 
this formalism as well, since the Born approximation fails 
when the range of the scattering potential R is too large 
in such a way that the associated energy scale hvpR^^ 
exceeds its strength, VR ^ hvp. The same behavior 
S ~ A/ei? is deduced (see Supplementary Information). 

Discussion and conclusions. As we have seen, the av- 
eraged amount of spin relaxed during a scattering event 
behaves as S* ~ A/e^?, independently of the nature of the 
scatterer, implying the general relation r^ « e|,Tp/A^. 
This result is not consistent with a linear scaling be- 
tween the spin relaxation time and the diffusion coeffi- 
cient at different gate voltages as it is observed in the 
experiments^ii^, suggesting that other mechanisms dom- 
inate spin scattering. This is consistent with the fact 
that CVD (Chemical Vapor Deposition) graphene-based 
spin valves show essentially the same spin transport prop- 
erties as exfoliated graphene^. This result implies that 
differences between exfoliated and CVD graphene, such 
as grain boundaries, do not limit spin transport. It is 



interesting to consider in detail the experimental data 
of Ref. \\M, The results show a sub-linear dependence 
of the diffusion constant on carrier density (proportional 
to the momentum scattering time) and also of the spin 
relaxation time. This is clearly inconsistent with our re- 
sult if one assumes the Elliot- Yafet induced by defects in 
graphene as the main spin relaxation mechanism. 

The main exception to the Elliot relation comes from 
impurities which enhance locally the SO coupling, such 
as heavy impuritie^i^, since in that case there are ad- 
ditional channels for spin relaxation. That is also the 
case of impurities that hybridize directly with graphene 
carbon atoms, such as hydrogeni^. In that case, an en- 
hancement of the SO coupling, Aioc(r) is induced due 
to the local distortion of the lattice coordination. For 
Ajoc ^ A it can be shown that a ex (Af^^) / Ef^^, where 
Eioc is an energy scale comparable to the local shift of the 
chemical potential in the region where the spin orbit cou- 
pling is modified. Alternatively, scattering by local spins 
can modify significantly the spin relaxation, whose effect 
could be determined by the dependence on temperature, 
magnetic field, or injected current. It is worth mention- 
ing that the combination of Zeeman coupling with local 
moments and enhanced spin-orbit coupling can lead to 
interesting new effects^^. 

Another interesting consequence of our work is that, if 
we suppose weak scatterers so that Tp oc ep , the Elliot- 
Yafet mechanism implies r^ oc ep, which scales with 
carrier density in the same way as t~^. This behavior 
makes it difficult to distinguish this mechanism and the 
D'yakonov-Pcrcl one (t^ cx t~^), which could explain the 
behavior of spin lifetime at high temperatures^. 

Acknowledgments: We appreciate useful discussions 
with R. Kawakami, B. J. van Wees, J. Fabian, V. I. Falko, 
B. Ozyilmaz, S. Das Sarma, and A. K. Geim. AHCN ac- 
knowledges DOE grant DE-FG02-08ER46512 and ONR 
grant MURI N00014-09-1-1063. HO acknowledges finan- 
cial support through grant JAE-Pre (CSIC, Spain). This 
work was also supported by MICINN (Spain) through 
grants FIS2008-00124 and CONSOLIDER CSD2007- 
00010. 



D. Huertas-Hernando, F. Guinea, and A. Brataas, Phys. 

Rev. B 74, 155426 (2006). 

N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and 

B. J. van Wees, Nature (London) 448, 571 (2007). 

W. Han and R. K. Kawakami, arXiv:1012.3435vl [cond- 

mat.mtrl-sci] (2010). 

T.-Y. Yang, J. Balakrishnan, F. Volmer, A. Avsar, 

M. Jaiswal, S. R. AH, A. Pachoud, M. Zeng, M. Popinciuc, 

G. Guntherodt, et al., arXiv:1012.1156vl [cond-mat.mes- 

hall] (2010). 

A. Avsar, T.-Y. Yang, S.-K. Bae, J. Balakrishnan, 

F. Volmer, M. Jaiswal, Z. Yi, S. R. Ali, G. Guntherodt, 



B.-H. Hong, et al., arXiv:1104.4715vl [cond-mat.mes-hall] 

(2011). 
^ M. I. D'yakonov and V. I. Perel', Sov. Phys. Solid State 13, 

3023 (1971). 
■^ P. G. Elliot, Phys. Rev. 96, 266 (1954). 
* Y. Yafet, Solid State Physics (Academic, New York, 1963). 
^ D. Huertas-Hernando, F. Guinea, and A. Brataas, Phys. 

Rev. Lett. 103, 146801 (2009). 
^° G. Ertler, S. Konschuh, M. Gmitra, and J. Fabian, Phys. 

Rev. B 80, 041405(R) (2009). 
" F. Guinea, New J. of Phys. 12, 083063 (2010). 
^2 J. N. Chazalviel, Phys. Rev. B 11, 1555 (1975). 



1^ p. H. Song and K. W. Kim, Phys. Rev. B 66, 035207 

(2002). 
" F. J. Jedema, M. S. Nijboer, A. T. Filip, and B. J. van 

Wees, Phys. Rev. B 67, 085319 (2003). 
^^ C. Jozsa, T. Maassen, M. Poponciuc, P. J. Zomer, 

A. Veligura, H. T. Jonkman, and B. J. van Weels, Phys. 



^® M. Hentschel and F. Guinea, Phys. Rev. B 76, 115407 
(2007). 



19 
20 



16 



Rev. B 80, 241403(R) (2009). 

A. H. Castro-Neto and F. Guinea, Phys. Rev. Lett. 103, 
026804 (2009). 

C. Weeks, J. Hu, J. Ahcea, M. Franz, and R. Wu, 
arXiv:1104.3282vl [cond-mat.mes-hall] (2011). 



D. S. Novikov, Phys. Rev. B 76, 245435 (2007). 
A. R. Akhmerov and C. W. Beenakker, Phys. Rev. B 77, 
085423 (2008). 
^^ M. I. Katsnelson, F. Guinea, and A. K. Geim, Phys. Rev. 

B 79, 195426 (2009). 
22 Z. Qiao, S. A. Yang, W. Feng, W.-K. Tse, J. Ding, Y. Yao, 
J. Wang, and Q. Niu, Phys. Rev. B 82, 161414 (2010). 



J 



Supplementary information 

A. Spin relaxation Avithin the Born approximation 

Let's consider the scattering problem defined by the Hamiltonian "Hq + U, where Hq is nothing but the free 
Hamiltonian of Eq. (4) in the main text, and C/ is a generic scattering potential, diagonal in sublattice and spin 
indices. We take A = for the moment. We treat the problem in perturbation theory. To first order in the scattering 
potential (Born approximation), the wave function solution of the problem with energy e = Hvpk reads "if^^' + '^^^' , 
where ^'^^^ is solution in the absence of the scattering potential, and ^'■"'^•' is given by: 



*(i) (r) = / d^r'G (e, r - r') \Ho + el] U (r') *<"' (r') 



(15) 



where 



G(e,r)= {e + iOf-Hl 



(16) 



In this case, the Green function (IT51) has a trivial structure in spin and sublattice indices, G (e, r) = G (e, r) I, where 
it reads: 



1 /■ piqr ■ 

{2T:f J {e + iOy -{hvFq) 4 (fivj.) 
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(17) 



which is nothing but the 2D Klein-Gordon propagator (iJ„ are the Hankcl functions of first kind). The asymptotic 
form of the Green function is: 



G(e,r)«- 
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In terms of a conventional scattering problem, we take as incoming wave function ^*" = "if^^^ a Bloch state with 
positive helicity. In the light of the asymptotic form of the Green function, it is clear that the scattered wave behaves 
as e^^"^ l\J —ir in the asymptotic limit. Moreover, as it is sketched in Fig. 1 of the main text, for each scattering angle 
0, two different scattered waves must be considered, with opposite helicity. Then, the wave function of the problem 
can be written in the asymptotic limit as follows: 
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where fh {0) is the scattering amplitude at each chiral channel. By assuming the usual approximation |r — r'| w r — r' -f , 
and integrating by parts is easy to see that: 
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where k' = fcf . Then it is straightforward to obtain the scattering amplitudes given in Eq. (6) of the main text. 



In order to compute exactly the scattering amplitudes in the presence of the SO interaction we must add the 
Rashba-like coupling to "Hq- The first consequence is that the Green function acquires a non-trivial structure in spin 
and sublattice degrees of freedom. After some tedious calculus we obtain: 

G (e, r) = -^77^ W'' (^+0 + H^o'^ (fc_r)) I + _ 'f^ ^^ (h^^^ {k+r) - H^'^ {k^r)] (l + a,® s. 
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(21) 



where k± = {hvp) \/ 1^ ^ eA, s = {sx ± JSj^) /2, and ct^ (s^) are Pauli matrixes associated to the sublattice (spin) 
degrees of freedom. We can repeat the same analysis as before. For an incoming wave with energy e in the positive 
quasi-helicity channel there are two outgoing radial waves with the same energy, with associated scattering amplitudes: 



/^ {B) = - {hvF\ 
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(22) 



where q± = k'_^ — k is the transferred momentum, with kj_ = k±r. Note that in the presence of the Rashba-like 
coupling the spin and sublattice degrees of freedoms are completely entangled. As consequence, backscattering is 
not forbidden, since eigenstates of Hq given by Eq. (5) of the main text are not longer eigenstates of the sublattice 
chirality operator ^. In fact: 

(o- X .s)^ , ■ 
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The results for weak and Coulomb scatterers were obtained by evaluating numerically Eqs. (7)-(8) of the main text 
with the expressions deduced here. 



B. Spin relaxation for scattering centers with cylindrical symmetry 

1. Strong scatterers 



As it is well known, scattering by strong scatterers such as vacancies or resonant impurities cannot be studied 
within the Born approximation. Strong scatterers are used to be described as infinite potentials with a range of the 
order of the lattice constant. We assume isotropy in order to exploit the cylindrical symmetry of the problem by 
using the decomposition of the eigenstates into partial waves with well-defined generalized total angular momentum 
J = Iz + <yz/'2 + Sz/2, which is actually a global symmetry of the problem, where Iz is third component of the orbital 
angular momentum operator Iz = —i{xdy — ydx)- For each incoming wave with energy e, there are two reflected 
waves with the same energy, characterized by the reflection coefficients r\, rj , where / labels the angular momentum 
channel J = I + 1. Then, following^, we define: 
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where r J 2 is nothing but r\ 2 in the absence of SO coupling. This quantity can be seen as the amount of spin relaxation 
during a scattering event. 

Neglecting mixing of the two inequivalent valleys, an incoming wave with energy e and total angular momentum 
J ^ I + 1 can be written as: 
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where c± 



> '^ , and k± is defined as above. The outgoing wave can be written as the superposition: 
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FIG. 4; Weak scatterers: S as a function of the carrier concentration for A = 1 meV (in blue) and A = 0.5 meV (in red), in 
both cases _R = 1 A and Vo = 0.1 eV. These results are obtained by solving numerically the system (|3H) . Inset: Dependence 
on fc^\ A = 0.5 meV. 

We model a strong scattcrer as a circular void of radius R. By imposing zig-zag boundary conditions at r=R, we 
obtain the following expressions for the reflection coefficients at each / channel: 



k-Ji+i (k+R) Yi (k-R) + k+Ji (k+R) Yi+i (k-R) 
k+Yi {k+R) y,+i {k_R) + k_Yi {k_R) Yi+^ {k+R) 
k^Ji+i {k+R) Yi {k+R) - k^Ji {k+R) Yi+^ {k+R) 
k+Yi {k+R) Yi+i {k^R) + k^Yi {k^R) Yi+i {k+R) 



(27) 
(28) 



At low energies kR <C 1, 5* is dominated by the first harmonics ^ = —1,0. In fact, the following asymptotic expression 
is deduced: 



5: 



■K A 

2 ' e^ 



(29) 



which fits very well the exact result, as it is shown in Fig. 3 of the main text. 



2. Clusters of impurities 



In the case of weak scatterers, if the range of the scattering potential R is too large in such a way that the associated 
energy scale hvpR~^ exceeds its strength, VR ^ Hvp, then the Born approximation fails. We employ the same model 
as in the case of strong scatterers. The incoming and outgoing waves are ((25|) and ([26]) respectively. Inside the 
potential r < R, the wave function regular at the origin can be written as the superposition: 
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ide 



Ji {q+r) e'"> 
tc'+Ji+,{q+r)e^('+')' 



®lt) 



+4 



„iW 



Ji {q-r) e' 
tc'_J;+i(g-r)e'('+i)'' 



®\t) 



c'+Ji+i{q+r)e'('+^^' 
iJi+2 {q+r) e'('+2)'' 

"c'_Ji+i(q-r)e'('+i)'' 
iJi+2{q-r)e^^'+^'>f> 



®|i) 



®u> 



(30) 



where q± = y^(e - y)^ =F A (e - F)/ {hvp) and c'+ = 



e-V 
hvFq± 



By imposing matching conditions at r = i? we obtain 4 equations that define the unknowns of the problem, in 
particular r\ and rj: 



c+Ji+i {k+R) 
c+Ji+i {k+R) 



Ji {k+R) 
r[c+Yi+i {k+R) 4 

r[c+Yi+i {k+R) - 



r[Yi {k+R) + r!,Yi {k^R) = t[Ji {q+R) + t'^Ji {q^R) 

' " ^- - ^'" '^^ - t[c'+Ji+i {q+R) + t'^c'_Ji+i {q-R) 

' "' ^- {q+R) 



r\cJYi+^ {k^R) 



ric^Yi+i (fc_i?) = t[c'+Ji 



+ ■^1+1 



t'2c'_Ji+i {q^R) 



Ji+2 {k+R) + r{Yi+2 {k+R) - r\Yi+2 {k^R) = t{Ji+2 {q+R) - tlJi+2 {q-R) 



(31) 
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FIG. 5: The reflection coefficient ri in the limit kpR ^ 1 for y = 0.1 eV and V — 0.5 eV, and R = 20 nm. ri tends to the 
asymptotic value — tan ( -^^ J , -0.033 and -0.167 respectively, as kp grows. 



In the absence of SO coupling, the reflection coefficients read {£, = sgn (e — V)): 



„M 



Ui (qR) Ji+i (kR) - Ji {kR) J,+i {qR) , ^Ji+, (qR) J1+2 {kR) - Ji+i (kR) J1+2 (qR) 



Yi {kR) Ji+i {qR) - 
^Ji{qR)Ji+i{kR) 



^Yi+i {kR) Ji {qR) 
- Ji {kR) Ji+i {qR) 



Yi+i {kR) Ji+2 {qR) - e^+2 {kR) Ji+i {qR)_ 

^Jl+l {qR) Jl+2 {kR) - Ji+i {kR) Ji+2 {qR) 

Yi {kR) Ji+i {qR) - ^Yi+i {kR) Ji {qR) Yi+i {kR) J1+2 {qR) - ^Yi+2 {kR) Ji+i {qR) 



(32) 



At low energies kpR <C 1, the most relevant channels correspond to / = 0,-1,-2. These reflection coefficients 



behave as 



VtR- 



^. In order to estimate the effect of the SO coupling in the doped regime it is enough to consider 



the energy splitting of the chiral sub-bands, as it is argued in the main text. The substitution e -> e ± A/2 -I- O (A) 



VeR^ 



AVR^ 



O (A/e)^ and ra 



VtR 



Y -|- O (A/e) . So the averaged amount of relaxed spin 



gives ri --.. ^^^^^2 , 2(rivFf ' ^^^i"' """" ' ^ '^ (Suf)"" 2(ftt), , 

behaves as S ^ ^/^i a-s it is numerically shown in Figure ID This is the same result {S is in fact of the same order) 

as the one deduced within the Born approximation. 

The opposite limit kpR 3> 1 has been studied in order to analyze scattering by clusters of impurities^. Assuming 
that charged impurities inside the cluster do not break the sublattice symmetry, so there is no gap opening, the main 
effect is a local shift of the chemical potential tp ^^ ^f + V (V > 0, note the change in the sign of V in relation with 
the previous calculation) inside the cluster. The effect of the long-range Coulomb potential has been analyzed before. 
The present approach is valid for circular shapes, where R is the radius of the cluster. For clusters large enough we 
can consider the regime characterized by kR ^1. If we take into account the asymptotic behavior of the Bessel 
functions: 



Ji{x) 



Yii^) 
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\ — cos X 7 
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it is easy to see that in this regime r2 ^ and: 



ri sa tan {Rk — Rq) 



fVR\ 



(33) 



(34) 



As it is showed in Fig. [Sj in this regime ri as a function of k presents some oscillations which attenuate as k grows, 
approaching the asymptotic value — ta.n {V R/ hvp) ■ These oscillations translate to the scattering cross section^, and, 
importantly, to S. However, these oscillations are expected to be averaged out for less symmetric cluster shapes, as 
it is deduced from semi-classical arguments^. In order to do so within our approach, we study the system (j3ip in the 
limit k±R^ 1. By taking the asymptotic form of the Bessel functions we arrive to: 



I c+ cos {q+R - f - ^) sin {k+R - ^ - l^) - c'+ cos {k+R - f - ^) sin {q+R - f - ^) 



c+ cos {q+R - 



\k+R 



cJj^ sin {q+R 



sin {k+R - 



l2L\ 
2 I 



(35) 
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FIG. 6: 5* as a function of A in the case of scattering by clusters of impurities for V — 0.5 eV, R = 20 nm and kp ~ 0.01 A ^ . 
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FIG. 7: S as a function of kp in the case of scattering by clusters of impurities for V — 0.5 eV, R = 20 nm and A = 0.5 meV. 



Then, we define the amount of spin relaxation as: 



S 



2^1 r 1 I 



E 



„0l\ 

/111 



(36) 



The results are shown in Figs. |6]and[7l The oscillations of S* as a function of kp persist, but clearly the enveloping 
curve behaves as r^ 1/kp. So 5 ~ A/e^, as the previous cases. 



C. Effect of a local enhancement of the SO coupling 



As it is argued in the main text, if the scattcrcr induces a local enhancement of the SO coupling, then the Elliot 
relation docs not hold since there are additional channels for spin relaxation apart from the one induced by momentum 
scattering. Since the source is the same, the presence of the scatterer, a correlation between the momentum relaxation 
time and the spin relaxation time is expected. A quantity a can be defined, like in the main text, as the spin- flip 
probability during a momentum scattering event. Nevertheless, two different mechanisms contribute to a in this case, 
the spin- flip induced by momentum relaxation (Elliot-Yafet), and the spin-flip induced by the local SO coupling. 
The latter manifests itself in the calculation within the Born approximation as a new contribution to the scattering 
amplitude in the positive helicity channel. If we assume a local enhancement of the SO coupling near the scattcrcr 
^loc |-j.^ (^g ^ ^^ g^g ^ perturbation to 'Ho, then we obtain this new contribution to scattering in the positive helicity 



10 



channel given by: 



ft"'-{e) = -{hv,)-'^^-^eA'-e-^'^ (37) 

where A^ "^ is the Fourier transformation of the local coupling evaluated at the transferred momentum q+ defined as 
before. 

If we assume that A'°^ ^ A, as it is argued in the main text, then we can neglect the Elliot- Yafet and study 
separately this new mechanism. We take A = for convenience, since then the spin up and spin down channels are 
asymptotically well defined. If we consider an incoming Bloch state in the spin up channel with energy e = hvpk, 
then the scattering amplitudes in the spin up and spin down channels read: 

h{e)^-ihvF)-\[^A'°He~^<> (38) 



Clearly a ex 4^ .If one assumes that the scattering center consists on a finite region where both a local shift of 

the chemical potential of the order of Eioc and an enhancement of the SO coupling of the order of A;oc are induced, 
then a ex Af^^/Ef^^, as it is pointed out in the main text. That is the case of heavy impurities. The analysis is more 
complicated in the case of resonant impurities which induces a local enhancement of the SO coupling due to the local 
distortion of the lattice coordination. The Born approximation fails in that case. By considering a model as the one 
described iii^ and including also a local shift of the chemical potential, then it can be shown that the scattering cross 
section has a non-monotonic behavior as a function of the carrier concentration. In other words, different regimes 
must be taken into account and the previous statement is not so solid. Another interesting point is that in the case 
of resonant impurities the spatial decay of the local SO coupling is governed by the dispersion relation of a bands^. 
To consider a step function for this coupling is probably a crude approximation which influences on the dependence 
of the spin relaxation time on the carrier concentration. 
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